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We present exact results on the partition function of the g-state Potts model on 
various families of graphs G in a generalized external magnetic field that favors or 
disfavors spin values in a subset Ig = {1, s} of the total set of possible spin val- 
ues, Z{G, q, s, V, w), where v and w are temperature- and field-dependent Boltzmann 
variables. We remark on differences in thermodynamic behavior between our model 
with a generalized external magnetic field and the Potts model with a conventional 
magnetic field that favors or disfavors a single spin value. Exact results are also 
given for the interesting special case of the zero-temperature Potts antiferromagnet, 
corresponding to a set-weighted chromatic polynomial Ph{G, q, s, w) that counts the 
number of colorings of the vertices of G subject to the condition that colors of ad- 
jacent vertices are different, with a weighting w that favors or disfavors colors in 
the interval Is- We derive powerful new upper and lower bounds on Z{G, q, s,v,w) 
for the ferromagnetic case in terms of zero-field Potts partition functions with cer- 
tain transformed arguments. We also prove general inequalities for Z(G,q, s,v,w) 
on different families of tree graphs. As part of our analysis, we elucidate how the 
field-dependent Potts partition function and weighted-set chromatic polynomial dis- 
tinguish, respectively, between Tutte-equivalent and chromatically equivalent pairs 
of graphs. 
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I. INTRODUCTION 



In this paper we continue our study of the g-state Potts model in a generahzed external 
magnetic field that favors or disfavors a certain subset of spin values in the interval Ig = 
{l,...,s}, on various families of graphs G We denote a graph G = {V,E) by its 

vertex set V and its edge ( = bond) set E. The numbers of vertices, edges, and connected 
components of G are denoted, respectively, by n{G) = n, e{G), and k{G). In thermal 
equilibrium at temperature T, the partition function for the Potts model on the graph G in 
this field is given hj Z = ^^{o- } ^^'^^ with the Hamiltonian 

(ij) p=l I 

where z, j, I label vertices of G, Oi are classical spin variables on these vertices, taking 
values in the set Iq = {1, q}, /3 = {ksT)'^, (ij) denote pairs of adjacent vertices, J is the 
spin-spin interaction constant, and 

( H for 1 < n < s 
Hp={ . (1.2) 

10 foTS + l<p<q 

Thus, for positive (negative) H, the Hamiltonian favors (disfavors) spin values in the interval 
Ig. This is a generalization of a conventional magnetic field, which would favor or disfavor 
one particular spin value. We denote as the orthogonal complement of Is in Jg, i.e., 
= {s + 1, q}, and we use the notation 

K = (3J , h = (3H , y = e^, v = y-l , w = e'' . (1.3) 

The physical ranges of v are v > for the Potts ferromagnet, and — 1 < < for the Potts 
antiferromagnet . 

It is very useful to have a general graph-theoretic formula for Z that does not make any 
explicit reference to the spins ai or the summation over spin configurations, but instead 
expresses this function as a sum of terms arisirig from the spanning subgraphs G' C G. This 



irisirig 

formula was derived and analyzed in Refs. j^, \4\ and is 



k(G') 

Z{G,q,s,v,w) = J2 ^'^""'^ n ""(^O (1-4) 



G'CG i=l 

where 



u^ = q- s + sw"" = q + s{w"^ -1) . (1.5) 
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This generalizes a spanning subgraph formula for Z in the case s = 1 due to F. Y. Wu [5|, l6[. 
In the special case H = 0, Eq. (I1.4p reduces to the cluster formula for the zero- field Potts 
model partition function 8|-[7|, denoted Z{G,q,v), namely 

Z{G,q,v)= ^ g'=(«') . (1.6) 

G'CG 

The original definition of the Potts model, (II. ip . requires q to be in the set of positive 
integers N+ and s to be a non-negative integer. These restrictions are removed by Eq. (11. 4p . 
Furthermore, Eq. (II. 4p shows that Z is a polynomial in the variables q, s, v, and w, hence 
our notation Z{G,q, s,v,w). If two graphs Gi and G2 are disjoint, then Z{Gi U G2) = 
Z{Gi)Z{G2) so, without loss of generality, we will usually restrict to connected G (although 
Eq. (II. 4p leads to consideration of disconnected spanning subgraphs G'). 

An important special case is the zero-temperature antiferromagnet, K = —00, i.e., v = 
— 1, and we denote 

Ph{G,q,s,w) = Z{G,q,s,-l,w) . (1.7) 

In this case the only contributions to Z are those such that no two adjacent spins have the 
same value. Thus, Ph{G, q, s, w) counts the number of proper g-colorings of the vertices 
of G with a vertex weighting that either disfavors (for < if < 1) or favors (for w > 
1) colors in the interval Ig- Here, a proper g-coloring is defined as an assignment of q 
colors to the vertices of a graph G subject to the condition that no two adjacent vertices 
have the same color. We have denoted these coloring problems as DFSCP and FSCP for 
disfavored or favored weighted-set graph vertex coloring problems 0]. The associated set- 
weighted chromatic polynomial constitutes a generalization of the conventional (unweighted) 
chromatic polynomial, which counts the number of proper g-colorings of a graph G. Recent 
reviews of chromatic polynomials include j9|-[ll|. 

As is evident from the discussion above, the model defined by Eq. (II. ip with Eq. (II. 2p is 
of interest both in the context of statistical mechanics and in the context of mathematical 
graph theory. It also has an application to certain frequency allocation problems in electrical 
engineering 



II. SOME BASIC PROPERTIES OF Z{G,q,s,w,v) AND Ph{G,q,s,w) 

In this section we discuss some basic results about Z{G, q, s, v, w) and Ph{G, q, s, w) that 
will be needed in our work. Applying the factorization 

m— 1 

w"'-l = {w-l)J2w^ (2-1) 
i=o 



3 



in Eq. (11. 4p with m = n{G[), one sees that the variable s enters in Z{G,q, s,v,w), and 
Ph{G, q, s, w) only in the combination 



t = s{w-l). (2.2) 

Since Ig ^ Iq, whence < s < g, and since w > for any real external field H, it follows 
that 

Um = q + s(w" - 1) > . (2.3) 

Therefore, for the ferromagnetic case v > 0, each term in the sum over spanning subgraphs 
in Eq. (II. 4p is nonnegative. For a given spanning subgraph G' C G, consisting of a sum of 
k{G') connected components G[, where i = 1, k{G'), the contribution to Z{G,q, s,v,w) 
in Eq. (II. 4p is the number of spanning subgraphs G' of a particular topology, Nq', times 
^e(G') Yl'^^ ■* which has the generic form 

fc(G') 

iVG,^;^(«') H . (2.4) 

Here 

fc(G') 

E ^(^^) = ^ • (2.5) 

i=l 

Since some of the components G'^ and G^- may have the same number of vertices, n(G-) = 
n{G'j), the product in Eq. (12.41) can also be written as Ylji'^vj)^^ , where rj takes on certain 
values in the set {1, ...,n} and the exponents pj are integers taking on certain values in the 
set {1, k{G')}. As a consequence of Eq. (12.50 . these satisfy the relation 

J2pjrj = n. (2.6) 

j 

Note that Um satisfies the identity 

Um{q,s,w) = w"'um{q,q- s,w~^) (2.7) 

where we have written function of its three arguments g, s, w. A given spanning 

subgraph G' corresponds to a partition of the total set of vertices depending on which edges 
are present and which are absent. The sum of the coefficients Nc of the various terms 
Ng' YljiurjY^ that multiply a given power v'^'^^'^ in Eq. (I1.4p is {^^q)-^) since this is the 
number of ways of choosing e(G') edges out of a total of e(G) edges. These satisfy the 
relation 

e[G) 



e(G')=0 
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This reflects the fact that there are 2^*^"-^^ spanning subgraphs of G, as follows from the 
property that these are classified by choosing whether each edge is present or absent, and 
there are 2^^^^ such choices. In mathematical graph theory, a loop is defined as an edge that 
connects a vertex to itself and a cycle is a closed circuit along the edges of G. In the following 
we restrict to loopless graphs. For any such n-vertex graph G, the terms in Z{G,q,s,v,w) 
proportional to v^, v^, and y^C^) can be given in general, as 

Z{G, q, s, V, w) = < + e{G)vu2u'l'^ + ... + v^^^^Un ■ (2.9) 

The partition function Z{G, q, s,v^w) satisfies the following identities 

Z{G, q, s, V, 1) = Z{G, q, 0, v, w) = Z{G, q, v) , (2.10) 

(where, as above, Z{G,q,v) is the zero-field Potts partition function), 

Z{G, q, s, V, w) = Z{G, q,q — s, v, w^^) , (2-11) 

(c.f. Eq. dMD) and 

Z{G,q,q,v,w) =w'' Z{G,q,v) . (2.12) 

Setting V = —1 in these identities yields the corresponding relations for Ph{G,q, s,w); for 
example, Eq. (12. lip yields 

Ph{G, g, s, w) = Ph{G, q,q-s, w'^) . (2.13) 

There are a number of equivalent ways of writing Z{G,q, s,v,w) as sums of powers of a 
given variable with coefficients depending on the rest of the variables in the set {q, s,v,w}. 
The basic spanning subgraph formula (11.41) is a sum of powers of v. A second convenient 
form in which to express Z{G, q, s, v, w) is as a sum of powers of w with coefficients, denoted 
as f3z,G,jiQy s,'^), which are polynomials in q, s, and v: 

n 

Z{G, q, s, v,w) = Y^ Pz,G,j{q, s, v) . (2.14) 
i=o 

The symmetry (12.111) implies the following relation among the coefficients: 

f3z,G,j{q, s, v) = Pz,G,n-M, q-s,v) for < J < n . (2.15) 
For the special case f = — 1, we write 

n 

Ph{G, q, s,w) = J2 Pg,M^ , (2.16) 
i=o 
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where 

PgM^^) = Pz,gM^s,-1) . (2.17) 

From f l2.15p . we have 

= /3G,n-,(g,g-s) forO<j<n. (2.18) 

We have proved further that Q 

(3z,G,n{Q^s^v) = ZiG,s,v) (2.19) 

and 

/3z,G,o(g,s,^) = ^(G,g-s,t;) , (2.20) 

so that for v = —1, f^coil, s) = P{G, q — s) and Pcnil, s) = P{G, s). Various general factor- 
ization results were also given in Ref. 0] for these coefficients l3z,G,jiQ, s,v) l3G,j{q,s), 
including the following: 

For 1 < j < n, Pz,G,j{<l) s,v) and (3G,j{(l, s) contain a factor of s . (2.21) 

For < j < n — 1, /3z,G,j{(l, s, v) and (3G,j{q, s) contain a factor {q — s) . (2.22) 

The minimum number of colors needed for a proper g-coloring of a graph G is the chromatic 
number, x{G)- A further factorization property is that 

x(G)-i 

/3G,n(?5 s) contains the factor Y\ ~ j) ^ (2.23) 

j=0 

and 

X(G)-1 

f^Gfiiq, s) contains the factor Y\ (q ~ s ~ j) ■ (2.24) 

3=0 

A third useful type of expression for Z{G, q, s, v, w) is 

n 

Z{G, q, s, v,w) = ^ az,G,n-jis, v, w) g""^ . (2.25) 
i=o 

With the notation 

aG,n-i(s, w) = az,G,n-jis, -1, w) , (2.26) 

we then have 

n 

Ph{G, q, s,w) = J2 «G,n-,(s, w) g"-^' . (2.27) 
i=o 
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This form is particularly convenient for comparisons with the conventional unweighted chro- 
matic polynomial P{G, q) = Ph{G, q, 0, w) = Ph{G, q, s, 1). 

For a graph G, the number of linearly independent cycles, c{G) (the cyclotomic number), 
satisfies the relation 

c{G) = e{G) + k{G) - n{G) . (2.28) 

A connected n- vertex graph with no cycles is a tree graph, T„, while a general graph with 
no cycles, which can be disconnected, is called a forest graph. We denote a graph G with no 
cycles as Gnc and define 

q'=^ , v' = - . (2.29) 

s s 

In Ref. 0] we proved that for such a cycle- free graph Gnc, 

Z{Gnc, g, s, V, w) = s^Z{Gnc, 1, v' , w) . (2.30) 

This relation allows us to obtain Z{Gnc,q, s,v,w) from Z{Gnc,qA,v,w) for any cycle-free 
graph Gnc- In particular, all of the results for Z{G, q, s, v, w) for various types of tree graphs 
calculated in Ref. jsl for s = 1 can be used to obtain the analogous results for general s. 

For a graph G, let us denote the graph obtained by deleting an edge e E E a.s G — e 
and the graph obtained by deleting this edge and identifying the two vertices that had been 
connected by it as G/e. The Potts model partition function satisfies the deletion-contraction 
relation (DCR) 

Z{G, q, v) = Z{G - e, g, v) + vZ{G/e, q, v) , (2.31) 
and, setting v = —1, the chromatic polynomial thus satisfies the DCR 

P{G, q) = P{G -e,q)- P{G/e, q) . (2.32) 



However, as we showed in Ref. j^, in general, neither Z{G, q, s, v, w) nor Ph(G, q, s, w) sat- 
isfies the respective deletion-contraction relation, i.e., in general, Z{G, q, s, w, v) is not equal 
to Z{G — e, q, s, w, v) + vZ{G/e, q, s, w, v). The only cases where this deletion-contraction 
relation holds are for the values s = 0, w = 1, and w = where Z{G, q, s, v, w) reduces to a 
zero-field Potts model partition function. In Section 



III. UPPER AND LOWER BOUNDS ON Z(G, q, s, v, w) FOR v>0 

In this section we derive powerful new two-sided upper and lower bounds for the gen- 
eralized field-dependent partition function of the ferromagnetic {v > 0) Potts model, 
Z{G,q,s,v,w) on an arbitrary graph G in terms of the zero-field Potts model partition 
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functions Z{G,ui,v) and Z{G,ui/w,v), where ui = q + s{w — 1) (c.f. Eq. (11. 5p ). These are 
especially useful because the zero-field Potts model partition function is considerably easier 
to calculate than Z{G,q, s,v,w). Throughout this section, it is understood that q > 0, 
< s < q, and v > 0. The former two conditions are obvious for our present analysis, while 
the latter will often be indicated explicitly. 

We first derive a lower bound for Z{G,q,s,v,w) for the range w > 1. To begin, we 
observe that, from its definition in Eq. (11.51) and factorization property (12.11) . Um satisfies 



m— 1 

Um = Q + s{w"^ — 1) = q + s{w — 1) 

j=0 

> q + s{w — 1) = Ui for w > 1 . (3.1) 

Substituting this inequality into the expression for Z{G,q, s,v,w) in Eq. (I1.4p in terms of 
contributions from spanning subgraphs G' C G, we have, for the same conditions 

Z{G,q,s,v,w) = J] ^^(a'j 

G'CG i=l 

> J2 ^'^^'^ {uif^''^ for v>0 and w>l . (3.2) 

G'CG 

But the expression on the second line of Eq. (13. 2 p is just the zero-field Potts model partition 
function given in Eq. (II. 6p with its argument q replaced by ui, namely Z{G,ui,v). Hence, 
we have derived a lower bound on Z{G, q, s, v, w): 

Z{G,q, s,v,w) > Z{G,ui,v) for t> > and w>l. (3.3) 

For the interval < w < 1, the inequality (13. ip is reversed: 

Um < ui for < < 1 , (3.4) 

and thus Eq. (13. 2p is replaced by 

Z{G,q,s,v,w) < ^ (mi)^'(^') for w > and 0<w<l. (3.5) 

G'CG 

Therefore, we obtain a second inequality, which is an upper bound: 

Z{G, q, s, V, w) < Z{G, Ui, v) for t> > and < w < 1 . (3.6) 

To derive two-sided inequalities, we make use of the symmetry relation (12. lip , which maps 
the interval w > 1 to the interval < w < 1 and vice versa. Let us start with the case 
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w > 1, for which we have proved the lower bound f l3.3p . Now, from the symmetry relation 
(12.1 ip we know that Z{G, q, s, v, w) = w"'Z{G, q, s, v, w) where s = q — s and w = w"^. Since 
w G [0, 1], we can apply our upper bound ( 13. 6p to Z{G, q, s, v, w), getting the inequality 

Z{G,q, s,v,w) < Z{G,ui,v) , (3.7) 

where 

Ul = q + s{w - 1) = q + {q - s)(w~^ - 1) = ^ . (3.8) 

Combining (13. 7p with (13. 3p . we derive the two-sided inequality 

Z(G,ui,v) < Z(G,q,s,v,w) <w'' Z(G,—,v) for > and ti; > 1 . (3.9) 

w 

For the interval < w < 1, by the same type of reasoning, we extend our upper bound (13. 6p 
to the two-sided inequahty 

w"" Z(G,—,v) < Z(G,q,s,v,w) < Z(G,Ui,v) for > and < w < 1 . (3.10) 
w 

As two-sided inequalities, these are powerful restrictions on the generalized field-dependent 
Potts model partition function in terms of zero-field Potts model partition functions with q 
replaced by Ui and Ui/w. 

We next prove some factorization properties of the upper and lower differences in these 
two-sided inequalities. First, if w = 1, then since Z{G,q, s,v,l) = Z{G,q,v) and Ui = q, 
it follows that the two-sided inequalities (I3.10p and (13. 9 p reduce to equalities, i.e., both the 
upper and lower differences vanish. Second, if f = 0, then the only contributions in the 
respective Eqs. (II. 4p and (II. 6p are from the spanning subgraph with no edges (called the 
null graph, A^„), so Z{G,q, s,0,w) = (ui)", and Z{G,q,0) = g", whence Z{G,Ui,0) = (ui)" 
and w"' Z{G,Ui/w,0) = (mi)". Hence, again, in this v = case, the inequalities (I3.10p and 
(13. 9p reduce to equalities and the upper and lower differences vanish. Third, if s = 0, then 
Z{G,q,0,v,w) = Z{G,q,v) and ui = q, so that Z{G,ui,v) = Z{G,q,v). Hence, if s = 0, 
then the lower difference in (13. 9p and the upper difference in (I3.10p vanish. Fourth, if w = 0, 
then Z{G, q,s,v,0) = Z{G, q — s,v) and Ui = q — s, so Z{G, Ui,v) = Z{G, q — s,v)] therefore, 
again, the lower difference in (13. 9 p and the upper difference in (I3.10p vanish. Together, these 
four results prove that the difference 

Z{G,q, s,v,w) — Z{G,ui,v) contains the factor w{w — l)sv . (3-11) 

Fifth, if s = g, then Z{G,q,q,v,w) = w"' Z{G,q,v) and Ui = qw, so w"^ Z{G,Ui/w,v) = 
w"" Z{G,q,v). Hence, if s = q, then the upper difference in (13. 9p and the lower difference in 
(I3.10p vanish. Combining this with the first two results above, we have shown that 

Z{G,—,v) — Z{G,q, s,v,w) contains the factor {w — l){q — s)v . (3.12) 
w 
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It is also useful to characterize the difference between the zero-field Potts model partition 
functions that constitute the upper and lower bounds in these two-sided inequalities fl3.9p 
and ( 13.10p . For an arbitrary graph G, we have 

w^Z{G,'^,v)-Z{G,m,v)= y M'^^'U^"^''^"'^'''^ - ll > (3-13) 

w ^-^ L J 

G'CG 

where G' is a spanning subgraph of G. Now the right-hand side of Eq. f l3.13p is nonzero 
only if G has at least one edge, and, in this case, the only nonvanishing contributions have 
n{G) - k{G') > 1. It follows that 

w Z{G,—,v)-Z{G ,ui,v) contains a factor vui(w - 1) . (3.14) 
w 

It is worthwhile to give some illustrations of these two-sided inequalities f l3.9p and fl3.10p . 
We first do this for tree graphs. For any n- vertex tree graph T„, ii w > 1, then the inequality 
reads 

ui{ui + v)^^^ < Z{Tn, q, s,v,w) < ui{ui + wv)^~^ for t> > and w>l. (3.15) 

where we have used Z(Tn, q, v) = q{q + t>)"~^. If w G [0, 1], then the inequality fl3.10p reads 

uiiui + wv}''-^ < Z{Tn,q,s,v,w) <ui{ui + v)"-^ for w > and < w < 1 . (3.16) 

(This example also shows how the apparent singularity at w = arising from the Ui/w 
argument in Z{G,ui/w,v) on the left-hand side of the inequality (13.100 is removed by the 
factor, yielding a nonsingular expression.) One gains further insight by calculating the 
differences between the polynomials that constitute the upper bound, the middle term, 
Z{Tn,q, s,w,v), and the lower bound for various tree graphs. For the path graph L2 and 
w > 1, the differences that enter in the two-sided inequality (13.151) are 

Mi(mi + wv) — Z{L2, q, s, v, w) = {w — l)(g — s)v > (3.17) 

and 

Z{L2, q, s, V, w) — ui{ui + v) = w{w — l)sv > . (3.18) 

For w G [0, 1] the differences that enter in (13.161) are obvious reversals of these, viz., ui{ui + 
v) — Z{L2, q, s, V, w) = w{l—w)sv > and Z{L2, q, s, v, w)—Ui{ui+wv) = {l—w){q—s)v > 0. 
For the path graph L3 and w > 1, the differences in (13.151) are 

mi(mi + wvY — Z{L3, q, s, f , w) = {w — l)(g — s)v 2ui + v{w + 1) > (3.19) 
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and 



Z{L^, q, s, V, w) — Ui{ui + vY = w{w — l)sv 2ui + v{w + 1) 



> 



(3.20) 



and similarly for w G [0, 1]. 

Among n-vertex tree graphs, the star graph S'„ has a particularly simple field-dependent 
Potts partition function, which was given in Ref. and is derived by a direct evaluation of 
the general formula (11.41) (for any v): 



n-l 



Z{Sn,q,s,v,w) = 

j=0 



n — 1 



j n—l—j 



(q - s) [ 



— s)\q + s{w — 1) + V 



+ SW 



q + s{w — 1) + wv 



n-l 



(3.21) 



Here j is the number of edges in a given spanning subgraph G', and the numerical prefactor 
("~^) in the first line of Eq. (I3.2ip is the number of ways of choosing j edges out of the total 
number of edges, n — 1, in Sn- For v > 0, substituting this result (13.211) into the two-sided 
inequahties (13.151) and (I3.16p . we can derive general formulas for the respective upper and 
lower differences. If tt; > 1 we find, for the lower difference in (13.151) . 



Z{Sn, g, S, f , w) - Ui{Ui + v) 



n-l 



SW I (Ml + Wm)" ^ - (mi + m)""^ 



n-l 
j=0 

sw{w — 
> . 



n — I 
j 



(mi)"-^~^ {w' -1) 



n-l 



n — 1 
3 



1=0 



(3.22) 



In the same way, if w G [0, 1], then the upper difference Ui{ui + v)"-^^ — Z{Sn, q, s, v, w) in 
(I3.16P is given by minus the right-hand side of Eq. (I3.22p . Similarly, if w > 1, then for the 
upper difference in (I3.15P we calculate 



Ui{ui + wv)'^ - Z{Sn,q, s,v,w) = (q-s) {ui + wv)"' - {ui + v) 



\n-l 



n-l 

{q- s){w -l)vY^ 
> . 



n — 1 
j 



iui) 



n-l-j yi- 



i=0 



(3.23) 
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Again, if w G [0, 1], then the lower difference Z{Sn,q, s,v,w) — ui{ui + wf)"~^ in (13.161) is 
given by minus the right-hand side of Eq. fl3.23p . 

For the circuit graph C„, if w > 1, the inequahty (I3.15P reads Z{Cn,ui,v) < 
Z{Cn, q, s, V, w) < Z{Cn, Ui/w, v). Using the fact that Z{Cn, q, v) = {q + v)"-' + {q — l)f 
we can write this exphcitly as 

(Ml + t;)" + {ui - l)i;" < Z(C„, g, s, v, w) < {ui + wv)" + (mi - U7)w"-^w" . (3.24) 

For 6*2 (which has a double edge), the differences that enter in this two-sided inequality are 

Z{C2, ui/w, v) - Z{C2, q, s, V, w) = {q- s){w - l)v{v + 2) > (3.25) 

and 

Z{C2, q, s, V, w) - Z{C2, ui, v) = w{w - l)sv{v + 2) > . (3.26) 

Similar illustrations of the general inequalities f l3.15p and fl3.16p can be given for L„ and Cn 
with higher values of n and for other families of graphs. 

For w > 1, we can prove a lower bound on Um that is stronger than (13. ip . To do this, 
we use the basic inequality that for real positive numbers a,, the arithmetic mean is greater 
than or equal to the geometric mean, i.e., 

-J]a,>[na.] " (3.27) 

(with equality only if aj = Ofc V j, k). Applying this to the sum Y7j'=o ^^at appears in 
the factorization relation Eq. (12. ip . we have, for all w > 0, 

m—l m— 1 l/m 

[n^l ™- (3-28) 

j=0 j=0 

Now rij^^ — ""^^5 where p = Y^^=i 3- Using the summation formula Y^^j=ij — n(n + l) 
we calculate that p = {m — l)m/2. Hence, the inequality (13.281) for w > can be written as 

m— 1 

J2 > mw^'^-^'>^^ . (3.29) 

j=0 

Since = q + s{w"^ — 1) = q + s{w — 1) Yl^=i ^ '^^ "^^^ ^se the lower bound fl3.29p 
to obtain a stronger lower bound on if w > 1 (but not if w G [0, 1), since in that case 
the prefactor {w — 1) is negative). Consequently, for w > 1, substituting (13.290 into the 
expression for Um, we derive the lower bound 

M^ = g + s(^/;'^-l) >g + ms(ty-l)M;(™-l)/2 for W > 1 . (3.30) 
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Clearly, this is an improvement over the lower bound (13. ip . Substituting this result into Eq. 
fll.4p with Um = Un(G'i)i "we thus obtain the following improved lower bound on Z{G, q, s, v, w) 
for w > 1 and the ferromagnetic range f > (where G' is a spanning subgraph of G): 



KG') ^ 

for w > 1 and f > . 

I J. \ I, / \ / 

G'CG 1=1 



Z{G,q,s,v,w)> Y,^"^^'^ n [q + HG',)s{w-l)w^''^^''^-^^/^ 



{3.31] 



Note, however, that in contrast with our previous lower bound (I3.2p . the right-hand side of 
this inequality cannot, in general, be written in terms of a zero- field Potts model partition 
function since the terms in the product depend explicitly on n{G'j). 



IV. SOME THERMODYNAMIC PROPERTIES 

The zero-field Potts model Hamiltonian Ti and partition function Z are invariant under 
the global transformation in which (Xj — )■ gai V i G V^, with g G Sq, where Sq is the 
symmetric (= permutation) group on q objects. In the presence of the generalized external 
field defined in Eq. (II. 2p . this symmetry group of "H and Z is reduced to the tensor product 

Sq-^Ss^Sq.s- (4.1) 

This simplifies to the conventional situation in which the external field H favors or disfavors 
only a single spin value if s = 1 or s = g — 1, in which case the right-hand side of Eq. (14. ip 
is Sq-i. For s in the interval 

2 < s < g - 2 , (4.2) 

the general model of Eqs. (II. ip and (II. 2p exhibits properties that are interestingly different 
from those of a g-state Potts model in a conventional magnetic field. For example, in the 
conventional case, at a given temperature T, if H » \ J\, the interaction with the external 
field dominates over the spin-spin interaction, and if h = f3H is sufficiently large, the spins 
tend to be frozen to the single favored value. In contrast, here, at a given temperature T, 
provided that s lies in the interval (14. 2p . if \H\ » \J\, this effectively reduces the model to 
(i) an s-state Potts model if if > 0, or (ii) a (g — s)-state Potts model if if < 0. In this 
limit, for given values of g and s and a given graph (say a regular lattice), there are thus, in 
general, four types of possible models, depending on both the sign of H and the sign of J. As 
an illustration of this, let us consider the case g = 5, s = 2 on (the thermodynamic limit of) 
a square lattice. For H = 0, the ferromagnetic version of the model has a first-order phase 
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transition, with spontaneous breaking of the 1S5 symmetry, at Kc = ln(l + y^) — 1-17, while 
the antiferromagnetic version has no finite-temperature phase transition and is disordered 
even at T = jg, 0]. For H > and H >> | J|, the theory reduces effectively to a two- 
state Potts model, i.e., an Ising model. Because the square lattice is bipartite, there is an 
elementary mapping that relates the ferromagnetic and antiferromagnetic versions of the 
model, and, as is well known, both have a second-order phase transition, with spontaneous 
symmetry breaking of the S2 ~ Z2 symmetry, at \Kc\ = \n{l + \/2) ~ 0.881 (where K = (5J), 
with thermal and magnetic critical exponents yt = 1, yh = 15/8, described by the rational 
conformal field theory (RCFT) with central charge c = 1/2. For H < and \H\ » \J\, the 
theory effectively reduces to a three-state Potts model. In the ferromagnetic case, J > 0, this 
has a well-understood second-order phase transition, with spontaneous symmetry breaking 
of the S3 symmetry, at Kc = ln(l + y/S) ^ 1.01, with thermal and critical exponents 
yt = 6/5, yh = 28/15, described by a RCFT with central charge c = 4/5 joi 0, 12|. In 
the antiferromagnetic case, J < 0, the model has no finite-temperature phase transition 
but is critical at T = (without frustration), with nonzero ground-state entropy per site 
S/kB = (3/2)ln(4/3) ~ 0.432 QE^. 

In particular, an interesting difference with respect to the g-state Potts model with a 
conventional external magnetic field appears in the case in which the spin-spin interaction is 
antiferromagnetic, i.e., J < 0. In the conventional case, there is competition between the two 
terms in the Hamiltonian, and resultant frustration. Here the situation is altered and depends 
on the chromatic number x{G) of the graph, li H > 0, then there is frustration if s < xiG), 
and this becomes increasingly severe as the temperature decreases, but if s > x{Cr), then 
this frustration is absent, because it is possible to satisfy the antiferromagnetic short-range 
ordering preferred by the spin-spin interaction while also satisfying the assignments of spin 
values preferred by the interaction of spins with the external field. (Of course, the presence 
of this field does have an effect in restricting the preferred range of values of the spins.) 
Similarly, ii H < 0, then there is frustration if (q — s) < x{G) but not if (q — s) > x{G). As an 
example, we may consider the case g = 5, s = 2 on (the thermodynamic limit of) a triangular 
lattice. For H > with H » \J\, the model reduces to an Ising model, and (i) if J > 0, this 
has a symmetry-breaking second-order phase transition at = (1/2) In 3 ~ 0.549, in the 
same universality class as on the square lattice, while (ii) if J < 0, there is frustration and, 
as a consequence, the model has no finite-temperature phase transition, but is critical at 
T = 0, with nonzero ground-state entropy S/ks — 0.323 [3]. For H < with \H\ >> |J|, 
the model reduces to a three-state Potts model, and (iii) if J > 0, this has a symmetry- 
breaking second-order phase transition at = ln[cos(27r/9) + y/S sin(27r/9)] 0.631 [15], 
in the same universality class as on the square lattice; while (iv) if J < 0, it has a weakly 
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first-order symmetry-breaking phase transition at Kc — —1.59 (id Il7l |. with a completely 
ordered ground state, reflecting the fact that the chromatic number of the triangular lattice 
is = 3. The more general case where \H\ is not >> \J\ encompasses a rich variety 

of thermodynamic behavior depending on the signs of H and J, the ratio of \H/J\, the 
values of q and s, the dimensionality of the lattice, and, in the antiferromagnetic case, the 
type of (/-dimensional lattice. Note that if J = 0, then (i) S/ks = Ins for if > 0, and (ii) 
S/ks = ln(g - s) for H < 0. 

Although a one- dimensional spin system (with short-ranged spin-spin interactions, as is 
the case here) does not exhibit any flnite-temperature phase transition, it can still serve 
as a worthwhile illustration of some thermodynamic properties. A simple example of this 
type is provided by our model on an inflnite one- dimensional lattice, with either free or 
periodic boundary conditions. We denote a reduced, dimensionless free energy per site as 
/ = lim„_^oo(l/''^) InZ. Then from our analysis above, we have, for the limits as n — ?■ oo of 
the line and circuit graphs, {L} and {C}, 

fi{L}, q, s, V, w) = f{{C}, q, s, v, w) = fmiq, s, v, w) = ln(Az,i,o,+) , (4.3) 

where Xz,i,o,+ is given below in Eq. (15. 8p . From this the various thermodynamic quantities 
such as the internal energy, speciflc heat, entropy, etc. can be calculated. Also, from this, one 
can obtain the function $({6*}, q, s, w) = lim„_j,oo Ph{Gn, q, s, w)^/" for the n — )• oo limits of 
Gn = Ln, Cn- The function $({(?}, g, s, w) generalizes the ground state degeneracy per site 
of the zero-temperature Potts antiferromagnet, W{{G},q) = lim„^oo -P(G', g)^/". We have 

$({L},g,s,w) = $({C},g,s,w) = <l>i^(g,s,w) = (Az,i,o,+)L.-i • (4-4) 

We note the following reductions of $i_d, which follow from the general identities given above: 

$i^(g,0,u;) = $iz)(g,s,l) = g- 1 , (4.5) 

$iz5(g,s,0) = g-s-l , (4.6) 

and 

<^iDiq,q,w) = w{q-l) . (4.7) 

With the ranges < s < g and w > understood,, d^io/dq > for the nontrivial interval 
g > 2, reflecting the greater freedom of color assignments with increasing g. Furthermore, 
d^in/dw > 0, as is clear from the original Hamiltonian formulation in Eqs. (11.11) and (11. 2p . 
The derivative d^in/ds > if w > 1, and d^iu/ds < if < w < 1, which follows from 
the fact that the external fleld favors (disfavors) spin values in if w > 1 (w G [0, 1)). Plots 
of as a function of g and w for flxed s are similar to the s = 1 results shown in Figs. 2-4 
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of Ref. [3|, except that the minimal value of q allowed is now s instead of 1, and the line for 
w = is now $({L}, g, s, 0) = g — s — 1 rather than g — 2. We proceed to give exact results 
for Z{G, q, s, v, w) and Ph{G, q, s, w) for several families of graphs. 



V. PATH GRAPH L„ 

The path graph L„ is the graph consisting of n vertices with each vertex connected to 
the next one by one edge. One may picture this graph as forming a line, and in the physics 
literature this is commonly called a line graph. We use the alternate term "path graph" 
here because in mathematical graph theory the line graph L{G) of a graph G refers to a 
different object (namely the graph obtained by an ismorphism in which one maps the edges 
of G to the vertices of L{G) and connects these resultant vertices by edges if the edges of 
G are connected to the same vertex of G). For n > 2, the chromatic number is x{Ln) = 2. 
In 0] we gave some illustrative calculations of Z{Ln, q, s,v,w). Here we present a general 
formula for this partition function. Let 

. q — s + V sw 



q — s w{s + v) 




(5.2) 
(5.3) 



s, = []) . (5.4) 



and 



Then 

Z{Ln, q, s, V, w) = ul Hifi (rz,i,o)""^ Si (5.5) 

and Ph{Ln,q, s,w) = Z{Ln,q, s, —l,w). It is straightforward to verify that our result for 
Z{Ln,q, s,v,w) satisfies the relation fl2.30p . We note that 

det(T2,i,o) = v{q + v)w , (5.6) 

independent of s, and 

Ti{Tz,ifl) = q - s + V + w{s + v) . (5.7) 
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The eigenvalues of T^,i,o are the same as the eigenvalues with coefficients of degree d = for 
the circuit graph C„ given in Eqs. (5.3) of Ref. 0], namely 



,o,± 



q — s + V + w{s + f ) ± 



Thus, we can also write 



{q — s + V + w{s + v)Y — 4:vw{q + v) 



1/2 



(5. 



Z{Cn, q, s, V, w) = Tr [(Tz,i,o)1 + (s - 1) + (q - s - l)v^ . 



(5.9) 



The graphs L„, C„, and, more generally, lattice strip graphs of some transverse width 
Ly and length = m are examples of recursive families of graphs, i.e., graphs Gm that 
have the property that Gm+i can be constructed by starting with Gm and adding a given 
graph H or, if necessary, cutting and gluing in H. For these graphs, Z{Gm, q, s, v, w) has the 
structure of a sum of coefficients that are independent of the length m multiplied by m'th 
powers of some algebraic functions. The results for transfer matrices for the case s = 1 in 
Ref. [2j] elucidated this structure for s = 1, and our calculation of Z{Cn, q, s, v, w) in Ref. |4| 
and Z{Ln, q, s, v, w) here elucidate this structure for general s. Note that, by Eq. fl2.9p . the 



term in Z(C„, q, s, v, w) of highest order in v is f = [q + s{w — l)Yl 



j=0 



Jf", part of 



which gives rise to the last two terms in Eq. (15. 9p . We note that for s = or ti? = 1, one can 
check that our expressions for Z{Ln,q, s,v,w) and Z{Cn,q, s,v,w) simplify, respectively, to 
Z{Ln, q, v) = q{q + v)^~^ and Z{Gn, q, v) = {q + f )" + {q — l)f Going from the case of 
sw{w — 1) = to sw{w — 1) 7^ 0, Z{Ln,q, s,v,w) expands from a sum of one power to a 
sum involving two powers, and Z{Gn, q, s, v, w) expands from a sum of two powers to a sum 
of four powers. 

As our exact solutions for Z{Ln, q, s, v, w) and Z{Gn, q, s, v, w) show, the field-dependent 
Potts partition functions Z{G, q, s, v, w) do not, in general, have any common factor. This 
contrasts with the case of the zero-field Potts partition function, which always has an overall 
factor of q. Similarly, in the f = — 1 special case defining the set-weighted chromatic poly- 
nomial, the resultant polynomials Ph{G,q,s,w) do not, in general, have a common factor. 
For special values of s, Ph{G, q, s, w) may reduce to a form with a common factor. The case 
s = (and the case w = 1) for which this reduces to the conventional chromatic polynomial 
is well-known; in this case P{G,q) has, as a common factor, Il^i'o^'^l? ~ j)- Similarly, for 
s = q, Ph{G, q, q, w) has this common factor multipled by w". For the special case s = 1 and 
for a connected graph G with at least one edge, it was shown in Ref. [sj] that Ph{G, q, 1, w) 
contains a factor (g — 1). However, it is not true that for a special case such 
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connected graph G with at least one edge contains a factor of (g — s). For example, using 
the elementary result 

Z{L2, q, s, V, w) = s{s + v)w'^ + 2s{q — s)w + {q — s){q — s + v) , (5.10) 

one sees that Ph{L2,q, l,w) = (g - l)(g - 2 + 2w), but Ph{L2,q,2,w) = 2w^ + A{q-2)w + 
{q — 2)(g — 3), which has no common factor. 



VI. COMPLETE GRAPHS 

The complete graph Kn is the graph with n vertices such that each vertex is connected 
to every other vertex by one edge. The chromatic number is xi^n) = n and the number of 
edges is e{Kn) = (g)- The (conventional, unweighted) chromatic polynomial is 

n-l 

P{K^,q) = l[{q-j) . (6.1) 

3=0 



For our later calculations, we will need our previous result for Ph(Kn, q, s, w) from Ref. 
which we mention here. We introduce a symbol xq = x6{x), where 6{x) is the step function 
from M —7- {0, 1} defined as 9{x) = 1 if x > and 6{x) = if x < 0. Our result is [4| 

n 

PhiKr,, g, s,w) = J2 I^kM^ s) w' (6.2) 

E=0 

where 

PK,A<l.s)=r][\{{s-i)\[ n (6.3) 

^ ^ i=0 i=0 

Here it is understood that if the upper index on either of the two products in Eq. (16. 3 p is 
negative, that product is absent, so that the first product is absent for £ = and the second 
one is absent for i = n. Note that 

/5i^„,£(g, s) = l3K„,n-e{q, q-s) , (6.4) 

in agreement with the general symmetry (12.181) . Substituting this in Eq. (16. 2p shows explic- 
itly that our result for Ph{Kn, q, s, w) satisfies the symmetry relation (12.130 . Note that 
is not a recursive family of graphs, so one does not expect Ph{Kn, q, s, w) to have the form 
of a sum of coefficients multiplied by powers of certain algebraic functions, and it does not, 
in contrast to Ph{Gn, q, s, w) for recursive famihes On such as Cn or L„. 
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The calculation of Ph{Kn, q, s, w) for the cases Ki and K2 = L2 are elementary. For 
K3 = C3 our general formula fl6.2p yields 

Ph{K3,q,s,w) = PiK3,s)w^ + 3sis-l){q-s)w^ + 3siq-s){q-s-l)w + P{Ks,q-s) (6.5) 

while for K4 we have 

Ph{K^, q, s, w) = P{K^, s)w'^ + As{s - - 2)(g - s)w^ + 6s{s - l){q - s){q - s - l)w'^ 
+ 4s(g - s){q - s - l){q - s - 2)w + P{IU, q - s) . (6.6) 

VII. p- WHEEL GRAPHS Wh^P^ = Kp + Cn-p 

The ]9-wheel graph Whn^ is defined as 

Wh';^^ = Kp + Cn-p, (7.1) 

i.e., the join of the complete graph Kp with the circuit graph Cn-p- (Given two graphs G and 
H, the join, denoted G + H, is defined as the graph obtained by joining each of the vertices 
of G to each of the vertices of H). (Here and below, no confusion should result from the use 
of the symbol H for a graph and H for the external field; the meaning will be clear from 
context.) The family of Whl^^ graphs is a recursive family. For p = 1, Whn^ is the wheel 
graph. The central vertex can be regarded as forming the axle of the wheel, while the n — 1 
vertices of the Cn-i and their edges form the outer rim of the wheel. This is well-defined for 
n > 3, and in this range the chromatic number is x(W^^n) = 3 if is odd and xi^^n) = 4 
if n is even. Although Kp is not defined for p = 0, we may formally define Whn'' = Gn- For 
the zero- field case, i.e., for the usual, unweighted chromatic polynomial and for an arbitrary 
graph G, 

P{Kp + G, g) = P{Kp, q)P{G, q - p) = qi^p)P{G, q - p) , (7.2) 
where q^m) is the falling factorial, defined as 

m—l 

3=0 

This result is a consequence of the fact that in assigning colors to the p vertices of Kp, one 
must use p different colors, and then, because of the join condition, one must select from the 
other q —p colors to color the vertices of G. In particular, for Wh^P\ this gives 

P{Wh^^\q) = PiKp,q)PiGn-p,q-p) 
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{q-l-pT-P + {q-l-p){-ir~P . (7.4) 



Note that, for arbitrary p, this chromatic polynomial consists of the prefactor times the sum 
of the (n — p)'th powers of iV^y^(p) ;^ = 2 terms. For p = 1, this number can be seen to be the 



Ly = 1 special case of a general formula in Eq. (3.2.15) of Ref. 18| for the join of Ki with 
a width-Ly cyclic strip. 

For the weighted-set chromatic polynomial, we generalize this coloring method as follows. 
Consider first Ki + G. There are two possible types of choices for the color to be assigned 
to the vertex of Ki. One type is to choose this color to lie in the set Is- There are s ways to 
make this choice, and each gets a weighting factor of w. For each choice, one then performs 
the proper coloring of the vertices of G with the remaining q — 1 colors, of which only s — 1 
can be used from the set Ig] this is determined by Ph{G, q — 1, s — l,w). The second type 
of coloring is to choose the color assigned to the Ki vertex to lie in the orthogonal set I^. 
There are {q — s) ways to make this choice, and since this is not the weighted set, there is 
no weighting factor of w. For each such choice, one then performs the proper coloring of the 
vertices of G with the remaining q — 1 colors, of which all s colors in the set Is are available, 
but only q — s — 1 colors in the orthogonal set I^ are available. This yields the result 

Ph{Ki + G, q, s, w) = swPh{G, q - 1, s - l,w) + {q - s)Ph{G, q-l,s,w) . (7.5) 

To calculate Ph{Kp + G, q, s, w) for a given graph G, one first carries out the proper coloring 
of Ki + G, using the result f l7.5p . One then joins the next vertex of Kp to Ki + G to get 
K2 + G, using the relation Ki + {K^ + G) = Kr+i + G and iteratively applies Eq. (17. 5p . One 
continues in this manner to carry out the proper coloring of the full join Kp + G. This yields 

p 

Ph{Kp + G, g, s,w) = Y, Z^kUI' ^) Pf^iG, q-p,s-i, w)w' . (7.6) 



Utilizing this coloring method, we have calculated Ph{Whn \ q, s, w) for arbitrary n. Let 
us define 

a{p, q, s, w) = q — s — {p + 1) + {s — l)w = q — {p + 1) + s{w — 1) — w (7.7) 

and 



^Wh(p)/,±{^^s,w) = ^ a{p,q,s-i,w)±[a{p,q,s-i,wf + Aw{q-p-l)]^^^ 



for < £ < p 

(7.8; 
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We note that for i = 0, these ^wh(p),e,±L3y ^-re equal to the v = —1 special case of \z,i,o,j 
given in Eq. (5.3) of our earlier Ref. [41] for the circuit graph with the replacement of q by 
q — p (and with j = 1,2 corresponding to ± here). This is in accord with the fact that the 
effect of the join of Kp with G is that the proper g-coloring of G can only use q — p oi the 
original q colors. We define two additional terms that do not depend on q or s, 

'^iyh(p),2p+3 = (7-9) 

and 

The total number of A's for Ph{Whn \ q, s, w) is thus 

A^PMH^h(p)),A = 2(p + 2) . (7.11) 

Note that in contrast to the unweighted chromatic polynomial of Whn \ where the number 
of A's, Np(^^rf^(p)yx = 2, is independent of p, here this number depends on p. In terms of these 
quantities, we find, for the weighted-set chromatic polynomial for Whn \ the result 



Ph{Wh^^\q,s,w) = J2PKUl^s) 



+ 



+ 



£=0 
r P 

XI I^KpM^ s){q-s-p + i- l)w^ 



w 



(7.12) 



This formula applies for integer p > I and also for p = if one sets /^Kp/iq, s) = for 
p = 0. It can be checked that for p = 0, Eq. (17.121) reduces to our result for Ph{Cn, q, s, w) 
given as the special v = —1 case of Eqs. (5.3)-(5.5) in Ref. [3]. It can also be verified 
that for p = 1 and s = 1, Eq. (17.121) reduces to the result given for this case in Eqs. 
(3.30)-(3.32) in Ref. |3']. Furthermore, since the graph Wh'l^ = Ki + K3 = K4, it follows 
that Ph{Wh^^\q,s,w) = Ph{K4,q,s,w). The symmetry (I2.13P is realized as follows: the 
summation on the first line of Eq. (17.121) goes into itself, while the sum of the expressions 
on the two subsequent lines of Eq. (17.121) transforms into itself with the replacement of w 
by in these expressions and the prefactor appearing overall. One could also study 
,q,s,v,w), but we have focused here on Ph{Wh''P\ q, s,w), since its calculation 
can be performed by combinatoric methods associated with the proper g-coloring condition. 
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We give some explicit examples of set-weighted chromatic polynomials Ph{Whn\q,s,w) 
obtained from our general formula ( I7.12P in the first appendix. 

Following our notation in Ref. and earlier works, the n — )■ oo limit of a family of 
n-vertex graphs Gn is denoted {G} and the continuous accumulation set of the zeros of 
Ph{Gn,q, s,w) in the complex q plane is denoted Bq. For recursive families of graphs, this 
locus is determined as the solution of the equality in magnitude of two (or more) A's of 
dominant magnitude, as a function of q (with other variables held fixed). The other loci By, 
etc. are defined in an analogous manner. These loci are typically comprised of curves and 
possible line segments. For studies of the n — oo limit of chromatic polynomials and their 
generalization to weighted-set chromatic polynomials, the locus Bg is of primary interest. 
Depending on the family of graphs, the locus Bq may or may not cross the real q axis. If it 
does cross the real q axis, we denote the maximum (finite) point at which it crosses this axis 
as qc- Extending our previous result for the p = case of {G} = {Wh^'''^ in Eq. (7.17) of 



Ref. B, we find the following result for general p: 

q^ = 2 + p + ^!^—^ f oi {G} = {Wh^^H andO <w <1 and 1 < s < p + 2 . 
1 + w 

(7.13) 

Regarding connections of this general formula to previously determined special cases, (i) for 
s = oi w = 1, this reduces to the result qc = 2 + p for the n — )■ oo limit of the chromatic 
polynomial P{Wh^^\ q) given in Eq. (22) of Ref. [lo']; (ii) for p = 0, this reduces to the 
result for the n — )• oo of Ph(Gn, q, s, w) given in Eq. (7.17) of Ref. j^, and (iii) for s = 1, this 
reduces to the result for the n — )■ oo limit of Ph(Wh^^\ g, l,w) given in Eq. (10.1) of Ref. [3| 
(with the obvious notation change {G} — > {Wh}). For the relevant interval < w < 1, the 
value of qc in Eq. fl7.13p is (a) greater than the value qc = 2 + p for the unweighted chromatic 
polynomial; (b) a monotonically increasing function of s for fixed w in this DFSCP interval; 
and (c) a monotonically decreasing function of w. These properties are consequences of the 
greater suppression of color values in the set Is as w decreases in the DFSCP interval, finally 
restricting the vertex coloring to use colors from the orthogonal set as w reaches 0. Thus, 
as w decreases from 1 to 0, gc increases continuously from 2+pto2+p + s. In contrast, 
the left-hand part of the boundary locus Bq changes discontinuously; as w decreases by an 
arbitrarily small amount below 1, the point on the left where Bg crosses the real q axis jumps 
discontinuously from q = p to q = p+s. This behavior is in agreement with the fact that 
in the two limits w = 1 and w = 0, Bg is comprised, respectively, of the unit circle centered 
at g = 1 + p and the unit circle centered at g = 1 + s + p. The change in the nature of the 
locus for s > 2 + p follows via the corresponding generalization of the analysis in Ref. Q to 
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VIII. EFFECT OF MULTIPLE EDGES IN A GRAPH 



Consider a loopless graph G = (y,E). Replace each edge with i edges joining the same 
pair of vertices and denote the resuhant graph as Gfe- Then the following is a theorem: 

Z{Gie, q, s,v,w) = Z{G, q, s,Vi,w), where f£ = (f + 1)^ — 1 . (8.1) 

Clearly, if f = 0, then Z{G, q, s,v,w) = {q — s + sw)", independent of the edge set E of G. 
Hence, in this case, the operation of replacing each edge by i copies of the edge has no effect 
on the partition function. This is seen at an analytic level via the property that if f = 0, then 
also ve = for any (positive integer) i. Further, for v = —1, where Z{G, q, s, v, w) reduces to 
the weighted-set chromatic polynomial Ph{G, q, s, w), the proper g-coloring constraint is the 
same regardless of whether a given edge is replicated or not, so again the replication does 
not affect this polynomial. In Eq. (18.11) . this follows because if f = — 1, then also V£ = —1 
for any (positive integer) i. Combining these results, we note that 

Vi — V contains the factor v{v + 1) . (8.2) 

Consequently, for any graph G with at least one edge (so that the operation of edge replication 
is not vacuous) and for positive integer i, 

Z{Gee,q, s,v,w) — Z{G,q, s,V£,w) contains the factor t>(t> + 1) . (8.3) 



IX. EFFECTS OF DELETION AND CONTRACTION OF EDGES 

As noted above in Sect. [TTl in general, neither Z{G,q, s,v,w) nor Ph{G,q, s,w) satisfies 
the respective deletion-contraction relation. It is of interest to investigate how these poly- 
nomials deviate frpm the deletion-contraction relation. A natural measure of this deviation 
for a graph G is 

[AZ{G,e,q,s,v,w)]DCR = Z{G,q, s,v,w)- Z{G-e,q,s,v,w)+vZ{G/e,q,s,v,w) . (9.1) 
We also define [APh{G, e, q, s, w)]dcr = [AZ{G, e, g, s, —1, w)]£)cr- We showed that Q] 

[AZ{G,e,q, s,v,w)]dcr contains the factor svw{w — 1) , (9.2) 
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and hence [/S.Ph{G,e,q,s,w)]DCR contains a factor of sw{w — 1). A particularly elegant 
general formula can be obtained for this deviation in the case of the family of star graphs, 
Sni i.e., graphs consisting of one central vertex with n — 1 other vertices, each of which is 
only connected to this central vertex. We find (for the nontrivial range n > 2) 

[AZ(^„, q, s, V, w)]dcr = svw{w - 1) [g + s{w - 1) + wv^"^ . (9.3) 



X. CYCLE MEASURE 

In view of the relation fl2.30p . one can define the following function, which serves as a 
measure of the presence of cycles in a graph G: 

[AZ{G, q, s, V, w)]cydes = Z{G, q, s, v, w) - s^'Z{G, q\ 1, f ', w) . (10.1) 

where q' and v' were defined in Eq. 02.291) . Clearly, this difference vanishes if s = 1, so, since 
it is a rational function, 

[AZ(G, g, s, w)]c2ycZes contains the factor (s — 1) (10.2) 

In Ref. 0] we derived the result 

TAv/r- M (s-1)m„i;" (s - l)(g- s + sw;")t;" 

[l\Z{Gn, g, s, V, w)\cydes = = . (10.3) 

s s 

This reflects the fact that C„ contains one cycle. 

Here we present another example of this difference function. Let us define a path graph 

with n vertices and each edge replaced by i edges joining the same adjacent vertices as L„ £. 

Note that L2,2 = G2- For L3_2 we calculate 

[AZ(L3_2, q, s, V, w)]cydes = s{2s^ + v'^s + Avs + v^)w^ 

+ 2s^{q - s)w'^ + 2s^{q - s)w + {q - s)(-2s^ + v'^s + 4vs + 2sq + v"^) . (10.4) 



XI. USE OF Z{G,q,s,v,w) AND Ph{G,q,s,w) TO DISTINGUISH BETWEEN 
TUTTE-EQUIVALENT AND CHROMATICALLY EQUIVALENT GRAPHS 



A. General 



Two graphs G and H are defined to be (i) chromatically equivalent if they have the same 
chromatic polynomial, and (ii) Tutte-equivalent if they have the same Tutte polynomial, or 
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equivalently, zero-field Potts model partition function. Here the Tutte polynomial T{G, x, y) 
of a graph G is defined as 

T(G,x,y)= 5^(x-l)'=(«')-^'(^)(y-ir(^') , (11.1) 

G'CG 

where G' is a spanning subgraph of G (and c(G') and k{G') were defined above as, respec- 
tively, the number of linearly independent cycles and the number of connected components 
of G'). This is equivalent to the zero- field Potts model partition function, via the relation 

Z(G, g, v) = {x- lf^^\y - 1)" T{G, x, y) , (11.2) 

where ?/ = f + 1 as in Eq. (11. 3p and x = 1 + {q/v). The Tutte polynomial is of considerable 
interest in mathematical graph theory, since it encodes much information about a graph. 
However, although it distinguishes between many graphs, there exist other pairs of graphs 
G and H that are different but have the same Tutte polynomial. An important property of 
our generalized field-dependent Potts model partition function Z{G,q, s,v,w) is that it can 
distinguish between many Tutte-equivalent graphs. Similarly, an important property of the 
weighted-set chromatic polynomial is that it can distinguish between many chromatically 
equivalent graphs. We study this further in this section. This property is true for all w and 
s values except the special values w = 1, w = 0, s = 0, and s = q, for which Z{G, q, s, v, w) 
is reducible to a zero-field Potts partition function (as well as the trivial case v = 0) and 
similarly for Ph(G, g, s, w). reducible to a chromatic polynomial. In Ref. 0] we proved that 
for any two Tutte-equivalent graphs G and H, 

Z{G,q, s,v,w) — Z{H,q, s,v,w) contains the factor s{q — s)vw{w — 1) . (11-3) 

In the following, we will generally phrase our analysis in terms of how the field-dependent 
Potts partition function distinguishes between Tutte-equivalent graphs; the special cases of 
the various expressions for v = —1 then show how the weighted-set chromatic polynomial 
distinguishes between different chromatically equivalent graphs. 

B. Tree Graphs 

A class of Tutte-equivalent (and, hence also chromatically equivalent) graphs of particular 
interest is comprised of tree graphs, generically denoted T„. For these, T(Tn,x,y) = 
so 

Z{Tn,q,v) = qiq + vr-' and P(T„, g) = g(g - 1)""^ (11.4) 
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Note that e(T„) = n — 1 (and a tree graph cannot have any multiple edges). There is only 
one tree graph with n = 1 vertex, one with n = 2 vertices, and one with n = 3 vertices. 
There are two different tree graphs with n = 4 vertices, namely the path graph, L4, and the 
star graph, >S'4. Enumerations of tree graphs with larger numbers of vertices are given, e.g., 
in Refs. 
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2l|. Let us consider two different n- vertex tree graphs (which thus have n > 4), 
denoted Gt and Ht. Since these have the same number of edges, inspection of the general 
Eq. (12. 9p shows that for the difference Z{Gt,q, s,v,w) — Z{Ht,q, s,v,w), not only the f° 
and terms, but also the v'^ terms cancel. Hence, 

Z(Gt, q, s, V, w) — Z{Ht, q, s, v, w) contains the factor f ^ . (11-5) 

We recall that Ig ^ Iq, so that < s < q, and that w > 0, as follows for any physical field 
H. These properties will be understood implicitly in the following. As preparation for the 
derivation of an inequality concerning Z{Gt,q, s,v,w) for Sn and L„ graphs, it is useful to 
give some explicit examples. Let us consider the two tree graphs with n = 4 vertices, namely 
S'4 and L4. In the following, we will usually omit the arguments q,s,v,w in Z{G,q, s,v,w) 
for brevity of notation. We have given exact expressions for Z(Sn) in Eq. (3.5) of Ref. 0] 
and for Z{Ln) in Eq. (15. 5 p above. For our present purposes, we focus on the expressions in 
terms of the spanning subgraph expansion. For S'4, this is 

Z^S^) = + 3vu2ul + Sv'^UsUi + v^Ui , (11-6) 

while for L4 we have 

Z{L4) = Ml + 3vu2ul + ^;^(2u3^^l + uj) + v^U4 . (11.7) 

The difference in the structure of the term proportional to f ^ arises from the differences in 
the spanning subgraphs with two edges in S'4 and L4. Hence, 

Z{S4} - Z{Li) = v^{u3Ui - ul) = v^s{q - s)w{w - if . (11.8) 

Since the last expression will appear as a factor in the differences Z{Gt) — Z{Ht) to be 
presented below, we give it a symbol: 

^ = s{q - s)v^w{w - if (11.9) 

and note that 

/i>0, (11.10) 

so that ZiSi) - Z{Li) > 0. 
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There are three different tree graphs with n = 5 vertices: 5*5, L5, and a graph that we 
denote as Y^, which has the form of a Y, with the vertical part made up of three vertices 
and two edges (shown in Fig. 1 of Ref. j^). The graph Yn is the generahzation of this 
graph in which the vertical part is comprised of n — 2 vertices forming a path graph Pn-2 (so 
that Y4 = 5*4). The spanning subgraph expansions for these graphs, in order of decreasing 
maximal vertex degree, are 

2(85) = ul + Avu2ul + Gv'^usul + Av'^u^ui + v'^u^ , 

(11.11) 

Z{Y5) =ul + Avu2ul + 2v'^{2u3ul + ului) + v'^^'iu^Ui + + w^s , (11-12) 

and 

Z(L^) = ul + 4:VU2ul + Sv'^^u^ul + ului) + 2t>^('U4'Ui + U3U2) + v'^u^ . (11.13) 
Thus, for the differences, we have 

Z{S^) - Z{Y^) = 2v'^{u3ul - ului) + v^{u4,ui - U3U2) 

= fi[2ui + v{w + 1)] , (11.14) 

^(5*5) — Z{L^) = 3v'^{u3ul — ului) + 2t>^(-U4-Ui — U3U2) 

= fx[3ui + 2v{w + l)] , (11.15) 

and 

^(Ys) - ^(Ls) = V'^{U3UI - ului) + V^{u4Ui - U3U2) 

= fi[ui + v{w + 1)] . (11.16) 

Now (remembering that < s < g and w > 0), for the ferromagnetic range v > 0, for 
nonnegative a and b, one has 

aui + bv{w + 1) >0 . (11-17) 

Hence, for the ferromagnetic case, each of the differences ^(5*5) — Z{Y^), Zi^S^) — Z^L^), and 
Z{Y^) — Z{L^) is non-negative. 
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From these explicit examples, one sees that the origin of these inequalities can be traced 
to inequalities among products of the m^'s. We proceed to prove two lemmas and then a 
general theorem. Our first lemma is 

Un-\Ui> Un-iUi for n > 2 and 2 < £ < n — 2 . (11.18) 

To verify this lemma, we expand and factor the given expression: 

Un-l Ul — Un-l Ui = s{q — s)w{l + w"'~'^ — W^~^ — 



s{q - 




-1) 




n-l-2 


1-2 


s{q - 








1=0 


i=o 



This lemma shows that the difference — that appears multiplying in Eqs. f lll.8[ 
f lll.l4p . f lll.lSp . and flll.l6p is nonnegative, and similarly that the difference M4M1 — u-iU2 
that appears multiplying in the last three of these equations is nonnegative. 

Differences of the form Z{Gt) — Z{Ht) for higher values of n involve differences of higher 
products of Ur factors, and there is an analogous inequality for these products. We prove this 
as a second lemma. Let us consider a generic term in Eq. (II. 4p . for the spanning subgraph 
G' = ©G^ with k{G') connected components, G[, each with n{G'j) vertices. This has the 
form (12. 4 p satisfying the relation (12. 5p . Our second lemma is, with £ = n — k{G') + 1, 

fe(G') 

ueu1~^ > Yl Unicr) for n>2 and 1 < £ < n, i.e., 1 < k{G') < n . (11.20) 
i=i 

For example, for the case n = 6, this lemma yields the inequalities U4UI > u^, U4u\ > U2, 
and U4UI > U4U2. This lemma is proved by the same method as Lemma 1. 

Combining the expression for Z{Sn,q,s,v,w) in the first line of Eq. (I3.2ip with our other 
results above, we have the following theorem: For the ferromagnetic case, 

Z{Sn,q, s,v,w) — Z(Tn,q, s,v,w) > for f > (11.21) 

for any tree graph T„. This is proved by applying the two lemmas above to the terms in the 
spanning subgraph expansions of these partition functions for Sn and a generic tree graph 
T„. In the second appendix we give further explicit results for differences of field-dependent 
partition functions for tree graphs with n = 6 vertices. 



28 



The difference Z{Y^) — Z{L^) in Eq. f lll.l6p (where we omit the arguments for brevity of 
notation) can also be understood using the recursive relation for n > 5: 

n-4 

n-4 

+v-\Y.w^)[ZiY,)-Z{L,)] , (11.22) 

j=0 

where Z{y^ — Z{L^ = Z^S^) — Z^L^) = ^ was given in Eq. flll.Sp . For the ferromagnetic 
range f > 0, each term on the right-hand side of Eq. (111.221) is nonnegative, and hence this 
proves the inequality 

Z{Yn,q, s,v,w) — Z{Ln,q, s,v,w) > for t> > . (11.23) 

Combining ffTLZTD and IHTT^ . we have 

Z{Sn, q, s, V, w) > Z{Yn, q, s, v, w) > Z{Ln, q, s, v, w) for > . (11.24) 



C. Properties of Graphs Intersecting in a Complete Graph 

One class of chromatically equivalent graphs consists of graphs whose chromatic polyno- 
mials can be shown to be equal by an application of the complete graph intersection theorem. 
We recall this theorem. Let us consider a graph G that has the property of being composed 
of the union of two subgraphs, G = G1UG2, such that Gi nG'2 = -^m for some m. In the rest 
of this subsection, we assume that G has this property. Then P{G, q) satisfies the relation 

This is sometimes called the complete-graph intersection theorem (KIT) for chro- 
matic polynomials. In contrast, in general, Ph{G,q, s,w) is not equal to 
Ph{Gi,q, s,w)Ph{G2,q, s,w)/Ph{Km,q, s,w). This equality holds only for the four val- 
ues w = 1, w = 0, s = 0, and s = q where Ph{G, q, s, w) reduces to a chromatic polynomial. 
As a measure of the deviation from equality, we define 

lATDUfr^ M - TDuin \ Ph{Gi,q,s,w)Ph{G2,q,s,w) n 1 ok^ 

[APh{G, g, s, w)]kit = Ph{G, g, s, w) — — r . (11.26) 

Ph{Km,q,s,w) 

Let us consider a graph with n = 5 vertices, denoted Glkl, comprised of a triangle 
K = with two line segments L, each of length one edge, emanating outward from two 
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vertices of the triangle. This graph LKL has n = 5, e = 5 and c = 1. A second graph, Gkll, 
also with n = 5, e = 5 and c = 1, is comprised of a triangle with a line segment two edges 
long emanating outward from one vertex of the triangle. These graphs are Tutte-equivalent, 
with A 

T{Glkl. X, y) = T{Gkll, x, y) = x\x + + y) , (11.27) 

or equivalently, 

Z{Glkl, q, v) = ZiGKLL, q, v) = q{q + v){q^ + + ^v^ + v^) . (11.28) 
It follows that these graphs are also chromatically equivalent, with chromatic polynomial 

P{Glkl. q) = P{Gkll. q) = q{q - l)\q - 2) . (11.29) 

In contrast, the field-dependent Potts partition function and the weighted-set chromatic 
polynomial successfully distinguish between these graphs. For the LKL graph we calculate 

Z{Glkl, q, s, V, w) = Z{Glkl, s, v)w^ + s{q - s){s + v){5s'^ + lOsv + 7v^ + 2v^)w'^ 



+ s{q — s) {q — s + 2)t>^ + (7g — s)f^ + 10s^(g — s — v) + 15sqv 



w 



+ s{q - s) (s + 2)v^ + (6g + s)v'^ + 10(g - s)^(s - v) + 15(g - s)qv 



w 



+ s{q - s){q - s + v) 5(g - sf + 10(g - s)v + Iv"^ + 2v^ 



w 



+ Z{GLKL,q - S,V) 



11.30) 



where Z{GLKL,q,v) = Z{GKLL,q,v) was given above in Eq. (111.281) . For the Gkll graph 
we calculate 

Z{Gkll, q, s, V, w) = Z{Gkll, s, v)w^ 



+ s{q - s) 



w 



+ s{q - s) - 10/ + 10s^(g - v) + sv{15q + 2v - v^) + 6qv^ + qV^ + + v'^ 
+ s{q - s) - 10(g - sf + 10(g - sf{q - v) + {q - s)v{15q + 2v- v"^) 



w 



+ 6qv'' + qv" + Av'^ + 



w 
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+ s{q - s) 5{q - sf + 15{q - sfv + 16(g - s)v'^ + 2{q - s)v^ + + v'^ 



w 



+ Z{GKLL,q- S,V) . 



;ii.3i) 



Thus, Z{GLKL,q,s,v,w) is not, in general, equal to Z{GKLL,q,s,v,w) and, taking the 
V — —1 special case, Ph{GLKL, Q, s, w) is not, in general, equal to Ph{GKLL, Q, s, w). For the 
differences, we find 



Z{Glkl, q, s, V, w) - Z{Gkll, q, s, V, w) = 



11.32) 



s{w - 1) + (1 + w)v'^ + 2wv + q + 2v 
and thus 

Ph{GLKL, q, s, w) - Ph{GKLL, q, s, w) = s{q - s)w{w - if s{w -l) + q- l- w 



(11.33) 

These calculations provide another illustration of how Z{G, q, s, v, w) can distinguish between 
different graphs that yield the same Tutte polynomial, and how Ph{G, q, s, w) can distinguish 
between different graphs that yield the same chromatic polynomial. 

In the context of graphs that can be decomposed into the union of subgraphs that intersect 
in a complete graph, it is also useful to calculate Z and Ph for the graph consisting of two 
Ks's meeting at a common vertex, ixi, denoted Gkk- This graph has n = 5, e — 6, and 
c — 2. The Tutte polynomial is T{Gkk, x, y) — {x -\- -\- yY, or equivalently 



so that 



Z{GKK,q,v) = —— ^^q{q + Sqv + 3v +v) 

Z{Ki,q,v) 



nOKK, q) = = q{q " lf{q - 



;il.34) 



;il.35) 



P{K^,q) 

Note that these polynomials factorize. This is not the case with Z{GKK,q,s,v,w) and 
Ph{GKK, q,s,w). We calculate 



Z{GKK,q,s,v,w) = Z{Gkk,s,v)w^ + s{q - s)(s + w)(5s^ + 13sv + 12v^ + 4v''^)w^ 



+ 2s{q -s) - + s^{5q - 6v) + sv{9q - v'^) + bqv'^ + {q + 3)v^ + v^^ 



31 



+ 2s{q -s) - 5{q - sf + {q - sf{5q - 6v) + {q - s)v{9q - v'^) + bqv^ + (? + 3)^;^ + i;^] w'^ 
+ s{q - s){q - s + v) b{q - + l?,{q - s)v + I2v'^ + Av^'^w 

+ Z{GKK.q-s,v) (11.36) 
and hence 

PKGkk. s, w) = P{Gkk. s)w'' + s{q - s){s - ifi^s - 8)w^ 
+ 2s{q - s){s- 1) bs{q - s) + s - 4g + 2 + 2s{q - s){q- s-1) bs{q - s) + {q - s) - Aq + 2 



w 



+ s{q-s){q-s-lf\b{q-s)-S^w + P{GKK,q-s) 

This example thus further illustrates how the factorization properties of the field-dependent 
Potts model partition function and set-weighted chromatic polynomial differ from those of 
the zero-field Potts model partition function and (unweighted) chromatic polynomial. 
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XII. CONCLUSIONS 

In this paper we have presented exact results on Z{G,q,s,v,w), the partition function 
of the Potts model in an external generalized magnetic field that favors or disfavors spin 
values in a subset Is of the full set Ig on various families of graphs G, and on Ph{G, q, s, w), 
the weighted-set chromatic polynomial. In particular, we have presented new general cal- 
culations of Z{G, q, s.v.w) for the case of path (line) graphs L„ and Ph{Whn\q, s,w) for 
p- wheel graphs Whn '' . Wc have discussed various features of our exact results for path, line, 
circuit, star, and complete graphs. We have derived powerful new upper and lower bounds 
on Z{G,q,s,v,w) in terms of zero-field Potts partition functions with certain transformed 
arguments. We have also proved inequalities for the field-dependent Potts partition function 
on different families of tree graphs. An important property of Z{G, q, s, v, w) is the fact that 
it can distinguish between Tutte-equivalent graphs, and similarly, Ph{G, q, s, w) can distin- 
guish between certain chromatically equivalent graphs. We have elucidated this property 
with our inequalities and with explicit calculations. Some new results for quantities such 
as f{{G},q,s,v,w), ^{{G},q, s,w), and qc defined in the n — )■ cxd limits of recursive graph 
famihes have also been given. 

There are a number of interesting directions for future study. One direction is to investi- 
gate additional graph-theoretic applications of Z{G, q, s, v, w) and Ph{G, q, s, w). Second, it 
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is clearly valuable to calculate Z{G, q, s, v, w) and Ph{G, q, s, w) for other families of graphs, 
in particular, recursive ones such as lattice strips and study their properties. We have de- 
scribed some of the interesting differences in thermodynamic behavior between our model 
with a generalized external magnetic field and the Potts model with a conventional magnetic 
field. It would be worthwhile to explore these differences further for various values of q, s, 
H, and temperature on lattices in > 2 dimensions, using methods such as series expansions 
and Monte Carlo simulations. 
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XIII. APPENDIX 1 

In this appendix we present some explicit expressions for the set-weighted chromatic 
polynomials Ph{Whn\q, s,w) obtained from our general formula (I7.12p . We begin with 
case p = 1, i.e., the family of wheel graphs, Whn"^ . The first nondegenerate case is n = 4, 
for which Wh^^'' = i^4, as noted (cf. Eq. fl6.6p ). For n = 5, we have 

Ph{Wh'i\q, s, w) = s{s - - 2){s'^ - 5s + 7)w^ + s{s - l)[5s'^ - 19(s - l)](g - 
+ 2s{s - l){q - s){5qs - 7q - + 3s + Q)w^ 
+ 2s{q - s){q-s- l)(-4g + 5qs - + 6 - 3s)w'^ 
+ s{q - s){q - s - l)[5(g - s)^ - 19(g - s - l)]w 

+ {q- s){q - s - l)(g -s-2)[{q- - 5(g -s) + 7] . (13.1) 
For n = 6, 

Ph{Wh''Q\q, s,w) = s{s - l){s - 2){s - 3)(s^ - 4s + 5)w^ 
+ 2s(s - l)(s - 2)[3s2 - ll(s - l)](g - s)w^ 

+ 5s(s - l)(g - s){3qs^ - 9qs + 7q- 3s^ + Ts^ - s - 5)w^ + 20s(s - lf{q - s){q - s - 1)W 
+ 5s(g - s)(g - s - l)(3g^s - 2q^ + Qq - 6qs^ - 5qs + 3s^ + 7s^ - 5 + s)w^ 
+ 2s(g - s)(g - s - l)(g - s - 2)[(3(g - s)^ - ll(g - s - l)]w 

+ (g - s)(g - s - l)(g - s - 2)(g - s - 3)[(g - s)^ - 4(g - s) + 5] . (13.2) 
For the family of p- wheel graphs with p = 2, i.e., Whn\ we first note that Wh\ = 

('21 

fsTs + ^2 = and Wh);' = K2 + K3 = K^, both of which are subsumed by our general 
formula (16.21) . For Wh^ our general result (17.121) yields 

Ph{Whf\q, s, w) = s(s - l)(s - 2)(s - 3){s'^ - 7s + 13)w^ 

+ 2s(s - l)(s - 2)(3s2 - 17s + 25)(g - s)m;^ 
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+ s{s - l)(g - s){15qs^ - 63qs + 67q - 15s'^ + 505^ - 12s - 59)w^ 

+ 4s(s - l)(g - s){q - s - l)[5s(g - s) - 8q + 13]w^ 

+ siq - s){q -s- l)(15g2s - ISg^ - 37qs + 55q - SOgs^ + 12s + ISs^ + SOs^ - 59)^^ 

+ 2s(g - s)(g - s - l)(g - s - 2)[3(g - s)^ - 17(g - s) + 25]w 

+ (g - s){q - s - l)(g - s - 2)(g - s - 3)[(g - s)' - 7(g - s) + 13] . (13.3) 

Explicit examples for higher values of p and n can be calculated in a similar manner from 
our general formula (17.121) . 



XIV. APPENDIX 2 



As discussed in the text, an important property of Z{G,q, s,v,w) is that it can dis- 
tinguish between different Tutte-equivalent graphs, and a similarly important property of 
Ph{G, q, s, w) is that it can distinguish between many chromatically equivalent graphs. Tree 
graphs provide a basic context in which to explore this property, since any two tree graphs 
are both chromatically equivalent and Tutte-equivalent. In the present text and in earlier 
work we have given Z{G, q, s, v, w) for tree graphs up to and including n = 5 vertices. Here 
we give results for tree graphs with n = 6 vertices. There are six such tree graphs, as shown 
in Fig. 6 of Ref. ^|, reproduced here as Fig. [T]for the reader's convenience, labelled (i) L„, 
(ii) Yq, (iii) isoYe, (iv) Hq, (v) Crg, and (vi) 5*6. These are listed in order of increasing max- 
imal vertex degree A. There are thus (2) = 15 differences of Z(G, g, s, f , w) polynomials for 
these six graphs. We list these below, as differences with respect to the graph with highest 
maximal vertex degree, 5*6 first, then with respect to the graph with next-highest maximal 
vertex degree, Gtq, and so forth. Furthermore, since the arguments are the same for all of the 
partition functions, we omit them, writing Z{G, q, s, v, w) — Z{H, q, s, v, w) = Z{G) — Z{H). 
The differences are 

Z{S^) - ZiCre) = /i [(3s^ + 3sv + v^)w^ 

+ {6s{q - s) + 3qv + v^}w + {3{q - + 3{q - s)v + v^} (14.1) 
Z{Se) - Z{Hq) = /i[(2s + v)w + {2(g - s) + v}]^ (14.2) 
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Z{Se) - Z{IsoYe) = (5s' + 6sv + 2v^)w^ 



+ {10s{q -s) + 5qv + 2v'^}w + {5{q - sf + 6{q - s)v + 2v'} 



(14.3) 



Z{Sg) - Z{Yq) = II (5s' + &SV + 2^2)^2 



+ {10s(g - s) + Qqv + 3w'}w + {5(g - sf + 6(g - s)v + 2f'} 



(14.4) 



^(S'e) - Z{Lq) = /X (6s' + Ss-u + 3v')w' 
+ {12s(g - s) + 7gv + 4v'}w + {o{q - sf + 8(g - s)v + 3v'} 

^(Cre) - ^(//e) = g + s(w - 1) + v g + s(w - 1) + -yw] 
= jjL s{s + v)uP' + {2s(g — s) + i'(q' + v)}w -\- {q — s){q — s -\- v) 



(14.5) 



(14.6) 



^(Cre) - Z{IsoYq) = (s + 'y)(2s + v)w'^ + (2s + 'y){2(g - s) + vjw 



+ (g-s + w){2(g-s) +w} 



(14.7) 



^(Cre) - ^(Fe) = 2{g + s(w - 1)} + v{w + l) q + s{w - 1) + v{w + 1) 



[s + v)(2s + v)w' + {4s(g - s) + 3gv + 2v'}w + (g - s + v){2(q - s) + 



(14.8) 



ZiCre) - Z{Le) = (s + 't;)(3s + 2v)w^ + {6s(g - s) + 4qv + 3v^}w 



+ (g-s + w){3(g-s) + 2u} 



(14.9) 



Z{Hg) - Z{IsoYq) = /i (s + v) W + {2s(g - s) + + (g - s + w)'] (14.10) 
Z(//6) - ^(le) = A* [(s + ^^)'«^' + {2s(g - s) + 2?^; + + (? - s + t;)'] (14.11) 



36 



Z{He) - Z{Le) = /i 2(s + f + {As{q - s) + 3qv + 2v'^}w + 2{q - s + 
Z{IsoYq) — Z(Yq) = fiwv{q + v) = s{q — + v)[w{w — 1)]^ 



Z{IsoYe) - Z{Le) = /i q + s{w - I) + v{w + 1) 



[s + v)w + q — s + V 



(14.12) 
(14.13) 
(14.14) 

(14.15) 



Z{Yq) - Z{Lq) = fi {s + v)W + {2s{q - s) + v{q + v)}w + {q - s + v)^ 

In addition to our theorems (lll.2ip and flll.23p . we observe that (for < s < g and w > 0) 
all of these differences are non- negative for the ferromagnetic range f > 0. 
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FIG. 1: Tree graphs with n = 6 vertices. 
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